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Constraints upon the CKM angle 2 from Belle and BaBar 
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The Belle and BaBar experiments have measured branching fractions and CP asymmetries in the charmless 
decay modes B° — > n + n~ , B° — > p tt^, and S° — > p + p~ ■ From these measurements, contraints upon the CKM 
angle 4> 2 can be obtained. These constraints consistently indicate that cf> 2 is around 100° . 



1. INTRODUCTION 

The Standard Model predicts CP violation to 
occur in B° meson decays owing to a complex 
phase in the 3x3 Cabibbo-Kobayashi-Maskawa 
(CKM) mixing matrix. This phase is illustrated 

ub ua 1 

= as vectors in the complex 
plane: the phase results in a triangle of nonzero 
height. One interior angle of the triangle, de- 
noted (f> 1 or j3, is determined from B° — > J/ijj K° 
decays. 1 Another interior angle, 4>i or a, is deter- 
mined from charmless decays such as B —> ir + ir~ , 
B° — > p + tt~ , and B° — > p + p~ . To determine (j) 2 
requires measuring time-dependent decay rates; 
here we present such measurements from the 
Belle [1] and BaBar [2] experiments. 

In neutral B meson decays, CP violation arises 
predominantly because of interference between a 
B° — > / decay amplitude and a B° — ► B — > / am- 
plitude. For the final states considered here, there 
are two decay amplitudes possible: a b^u "tree" 
and a b— >d "penguin" (see Fig. 1). Because these 
amplitudes have different weak phases, additional 
information is needed to determine <p 2 , such as 
the size of the penguin amplitude or the differ- 
ence in strong phases between the penguin and 
tree amplitudes. 



by plotting the unitarity condition v ub v ud 

v; b v cd + v t * b v td 





Figure 1. Tree-level diagram (top) and penguin 
diagram (bottom) for B° — > tt + tt~ , B° — > p + ir~ , 
and B°^p + p~ decays. 



1 Charge-conjugate modes are included throughout this 
paper unless noted otherwise. 
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2. ANALYSIS 

The analyses of B° tt+tt--, B° 
and £?° — > decays have several similarities. 

Events are selected by requiring two opposite- 
charge pion-candidate tracks originating from the 
interaction region, and appending zero, one, or 
two 7r°'s. The charged pion identification criteria 
are based on information from cither a DIRC de- 
tector (BaBar) [3] or time-of-flight counters and 
aerogel cherenkov counters (Belle) [4]. Both ex- 
periments also use dE/dx information from the 
central tracking chamber. 

B decays are identified via two kinematic vari- 
ables: the "beam-constrained" mass, m bc , and 
the energy difference, AE. The former is defined 
as 'El — p 2 B and the latter as E B — E bl where p B 
is the reconstructed B momentum, E B is the re- 
constructed B energy, and E b is the beam energy, 
all evaluated in the e + e~ center-of-mass (CM) 
frame. After selection cuts, the m bc and AE dis- 
tributions are jointly fit for the signal event yields. 
This fit includes contributions from backgrounds, 
whose m hc -AE distributions are obtained from 
either Monte Carlo (MC) simulation or extrapo- 
lation from m hc -AE sidebands. 

A tagging algorithm is used to identify the fla- 
vor of the B signal decay, i.e., whether it is B° 
or 6°. This algorithm examines tracks not asso- 
ciated with the signal decay to identify the flavor 
of the non-signal B. It depends predominantly 
on identifying leptons or kaons. The signal-side 
tracks are fit for a signal decay vertex, and the 
tag-side tracks are fit for a tag-side decay vertex; 
the distance Az between vertices is to good ap- 
proximation proportional to the time difference 
between the B decays: Az w (/3 / yc)At, where 
/?7 is the Lorentz boost of the e+e~ system and 
equals 0.43 (0.56) for Belle (BaBar). One sub- 
sequently does an unbinned maximum likelihood 
(ML) fit to At to measure or constrain <j> 2 . 

The dominant background for all three de- 
cays is e + e~ — > qq continuum events, where 
q = u,d, s, c. To distinguish such events from 
e + e~ — > BB events, the event topology is used: in 
the CM frame, continuum events tend to be colli- 
mated along the beam directions while BB events 
tend to be spherical. In Belle, the "shape" of an 



event is typically quantified via Fox- Wolfram mo- 
ments [5] of the form h e — J2i. 3 PiPj Pe{ cos @ij)i 
where i runs over all tracks on the tagging side 
and j runs over all tracks on either the tagging 
side or the signal side. The function P e is the 
£th Legendre polynomial and 8^ is the angle be- 
tween momenta p i and p*. These moments are 
combined into a Fisher discriminant [6], and the 
discriminant is subsequently combined with the 
probability density function (pdf) for the cosine 
of the angle between the B direction and the elec- 
tron beam direction. This yields an overall likeli- 
hood C, which is evaluated for both a BB hypoth- 
esis and a continuum hypothesis. Signal B — > / 
events are separated from continuum events by 
cutting on the likelihood ratio L BB - / \C BB + Cqq) ■ 
In BaBar, B — ► / signal is separated from con- 
tinuum background using several methods. For 
B° — ► 7r + 7r~, a cut |cos6> sph | < 0.8 is imposed, 
where 8 sph is the angle between the sphericity 
axis of the B candidate and that of the rest of 
the event. A Fisher discriminant (J-) is then con- 
structed from J2iPi an d ^2iPi\ cos #i| 2 ) where p i 
is the momentum of particle i, ^ is the angle be- 
tween p i and the B thrust axis (both evaluated in 
the e + e~ CM frame), and i runs over all particles 
not associated with the B decay. A pdf for T is 
included in the ML fit to At. For B° -> p+ir^ 
and one [7] of two B° — > p + p~ analyses, a neu- 
ral network is used that includes the two event- 
shape variables from T . The output of the neural 
network is included in the At fit. For the other 
B° — > p + p~ analysis [8], a cut |cos# th | < 0.8 is 
made, where th is the angle between the thrust 
axis of the B candidate and that of the rest of the 
event. The analysis subsequently uses a Fisher 
discriminant constructed from II obscrvables. 



3. B°^7V + 7V- 

The decay time dependence of B° j B° -^ir + Tr~ 
decays is given by [9] 



dN 
dAt 



,-At/T 



l-qC n7r cos(AmAt) 

+qS 7r7T sin(AmAt) 



. (1) 
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where q — +1 (q = — 1) corresponds to B° (B°) 
tags, and Am is the B°-B° mass difference. The 
parameters and S V7r are CP- violating and re- 
lated to 4> 2 via [10] 
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i sin (5 
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(4> 1 + (f) 2 ) cos 5 - 



(3) 



(4) 



where X = (23.2±}j)° [11], |P/T| is the mag- 
nitude of the penguin amplitude relative to that 
of the tree amplitude, and S is the strong phase 
difference between the two amplitudes. If there 
were no penguin contribution, P = 0, C 7rn = 0, 
and S njr = sin 20 2 . Since Eqs. (2) and (3) have 
three unknown parameters, measuring and 
iS X7r determines a volume in (fr 2 - 5 - \P/T\ space. 

The most recent Belle measurement of C W7r and 
is with 140 ftr 1 of data [12]. Candidates 
must satisfy 5.271 GeV/c 2 < m bc < 5.287 GeV/c 2 
and AE< 0.064 GeV; the final event sample con- 
sists of 224 B°^TT+n- candidates and 149 B° — > 
7r + 7r~ candidates after background subtraction. 
The ratio of signal to background is ~ 0.3. These 
events are subjected to an unbinncd ML fit to At, 
in which additional pdf 's and resolution functions 
are included to account for backgrounds. There 
are two free parameters in the fit, and the results 
are = -0.58 ± 0.15 (stat) ± 0.07 (syst) and 
= -1.00 ± 0.21 (stat) ± 0.07 (syst). These 
values are consistent with previous Belle measure- 
ments [13] and indicate large CP violation. Fig. 2 
shows the At distributions for the q = ±1 sam- 
ples; a clear difference is seen between the distri- 
butions. Many cross-checks have been done for 
this analysis, including an independent "blind" 
analysis using a binned ML fit. The latter results 
are very close to those of the main fit. 

The Belle values for and prescribe a 
95% C.L. volume in (j> 2 - S - \P/T\ space. Slicing 
this volume at fixed \P/T\ gives a 95% C.L. con- 
straint in the 4> 2 -5 plane; slicing this volume at 




Figure 2. Belle results for B° — > 7r + n~ [12]: the 
At distributions of q = 1 tags (top), q = — 1 
tags (middle), and the resulting CP asymmetry 
(bottom). The smooth curves are projections of 
the unbinned ML fit. 



fixed 5 gives a constraint in the <fi 2 -\P/T\ plane. 
Two such projections are shown in Figs. 3a and 
3b; the resulting constraints are 90° < <f) 2 < 146° 
for \P/T\ < 0.45 (as predicted by QCD fac- 
torization [14] and perturbative QCD [15]), and 
\P/T\ > 0.17 for any value of S. 

The BaBar experiment has also measured C vv 
and using an unbinned ML fit [16]. The most 
recent result is from 205 fb _1 of data [17]; the val- 
ues obtained are C n7r = —0.09 ± 0.15 (stat) ± 
0.04 (syst) and = -0.30 ± 0.17 (stat) ± 

0.03 (syst). These values are inconsistent with the 
Belle result at the level of 3.2er [11]. The BaBar 
analysis differs from that of Belle in that fewer 
cuts are made to enrich the data sample; rather, 
additional pdf's for the discriminating variables 
are included in the likelihood function. A total of 
68 030 events are fit, and a signal yield of 467 ±33 
B° — >ir + ir~ decays is obtained. There are 46 free 
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Figure 3. Belle results for B° -^tt + it~: (a) con- 
straints in the 4> 2 -S plane for \P/T\ < 0.45; and 
(b) constraints in the (f> 2 - P/(\P\ + \T\) plane for 
all values of S. The dark blue region corresponds 
to la C.L., the light blue region to 90% C.L., and 
the green region to 95% C.L. 



parameters (including and S nv ) in the fit. 

The Belle and BaBar results can be averaged 
together to constrain </> 2 . However, such a con- 
straint requires knowledge or assumptions about 
\P/T\ or 5. A model based on SU(3) sym- 
metry (and including an 5P(3)-breaking factor 
fx/fir f° r tree amplitudes) and the measured 
rates for B° -> K°tt + and B° -> K+tt' indi- 
cates (f> 2 = (103 ± 17)° [18] (note: this uses the 
BaBar result for 113 fb" 1 of data). A preferred 
method is to use isopin symmetry and the mea- 



sured rates for B + — > 7r + 7r°, B° 



7r°7r°, and 



charge-conjugates; this method can determine <p 2 
with little theoretical uncertainty [19]. However, 
the decay B° — » ir°ir° has only recently been 
observed and the asymmetry between B° and 
B a measured [20]. Future measurements with 
higher statistics should yield an interesting con- 
straint on 2 . The overall (B° + B°) ->■ tt°tt 
branching fraction can be used to obtain an up- 
per bound [21] on the angular difference 9 = 
<t>2 -02 off. where s ™ = sin20 2eff (i.e., 2cff ->0 2 
as P— >0). Using C n7r and the most recent values 
of the above branching fractions [11] fluctuated 
by la in the conservative direction, one obtains 
9 < 36°. 



4. B l 



<■ /9+7T 



For B° — > p + ir~, the final state is not a CP 
eigenstate. There are thus four separate decays 
to consider: B° — * p ± 7r T and B° ^ p^ir^ . The 
decay rates can be parametrized as [22] 

cx (l±Ag P ) x 



dN(B^p ± n- 



dAt 



l-q(C piT ±AC piT )co S (AmAt) 



q(S p7T ±AS p7T )sm(AmAt) 



(5) 



where q = +1 (g = —1) corresponds to B° (B°) 
tags. The parameters C p7T and S p7T are CP- 
violating, while the parameters AC plT and A<S 
are CP-conserving. AC p7r characterizes the differ- 
ence in rates between U W — > p" processes P° — » 
p + 7r~ or B° ^ p~n + and "spectator — » p" pro- 
cesses P° — > p~tt + or P° — > p + 7r~ (see Fig. 1). 

in addition, on differences in 



A<S p7r depends 
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phases between W — > p and spectator — > p am- 
plitudes. 

The parameter ^4^> is equal to the time and 
flavor integrated asymmetry: T(B — > p + 7r~) + 



r(B°- 



> p + 7T~ 



-) -r(B° -» p-7r+) -r(B° -> p- 



divided by the sum of the four rates. We 
define two separate CP asymmetries: 



7T+) 

also 



iV(fl°- 



>p-ir+)+N(B°^p+Tr-) 
C„„ + A p * ■ AC n 



1 + AC. 



/97T 



4 p7r • f7 



(6) 



and 



N(B°^p+tt-) - N(B°^p-7r+) 
N(B°^p+7T-) + N{B°^p-ir+) 



A p7r 

A CP 



C 



A P CP ■ AC p7T 



1 - AC, 



f)7V 



A p7r ■ r 

■"■CP ^pn 



(7) 



A^ depends only on W^p processes and A h 

depends only on spectator processes. 

Both BaBar and Belle have done unbinncd ML 
fits to the At distributions of .B — >p ± 7r ± decays 
to determine A P J P , C p7r , S p7r , AC p7r , AS p7I , A + _, 

and A , . The Belle analysis is with 140 fb _1 

of data [23]; the BaBar analysis, originally with 
81 fb _1 of data [24], has been updated with 
113 fb- 1 [25]. 

To remove charge-ambiguous decays and pos- 



sible interference between B° 



• /9+7T 



and B° ■ 



p 7r + amplitudes, one must eliminate the overlap 
region of the 7r + 7r~7r° Dalitz plot. Belle does this 



by requiring 0.57 GeV/c 2 < 
and to^o > 1.22 GeV/c 2 



T o<0.97 GeV/c 2 
BaBar makes 



the looser selection 0.40 GeV/c 2 < m i±]t0 < 
1.30 GeV/c 2 and requires that m iTi0 not be in 
this range. In addition, BaBar requires that the 
bachelor track from B — > pir has p > 2.4 GeV/c, 
where p is evaluated in the e + e~ CM frame; only 
14% of pions from (selected) p ± decays satisfy this 
requirement. Finally, Belle requires m w+w _ > 
0.97 GeV/c 2 to avoid the overlap region between 
B°^p°7r° and B° -> p ± n^ . 

Belle subsequently defines a signal region 
m bc > 5.27 GeV/c 2 ' and -0.10 GeV < AE < 
0.08 GeV. There are 1215 events in this region 
that pass all selection requirements. Fitting to 
the m hc -AE distributions yields 329 B° —> p + n~ 




Figure 4. Belle results for B° — > p + ir~ [23]: the 
At distributions of q = 1 tags (left), q = — 1 tags 
(right), and the resulting CP asymmetry (bottom 
rows). The asymmetry is shown separately for 
high-quality tags (r > 0.5) and low quality tags 
(r < 0.5). The smooth curves are projections of 
the unbinncd ML fit. 



candidates. The resulting At distributions for 
q = ±1 tagged events are shown in Fig. 4 along 
with projections of the unbinned ML fit in At. 
Also shown is the CP asymmetry which is con- 
sistent with zero. 

The BaBar results are similar to those from 
Belle; the corresponding At distributions and CP 
asymmetry are shown in Fig. 5. All Belle and 
BaBar results are listed in Table 1. There is very 
good agreement between the measurements ex- 
cept for AS pn , where the disagreement is ~ 2a. 
A recent BaBar analysis with 192 fb _1 of data [26] 
uses a different strategy than the quasi-two-body 
approach: it takes advantage of interference in 
the 7r + 7T~7r° Dalitz plot as prescribed in Ref. [27]. 
These results are also listed in Table 1 for com- 
parison; they are very similar to those from the 
quasi-two-body analyses. 

These measured values can be used to constrain 
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Table 1 



Results of fits to the At distributions for B°^p + n candidates. 





Belle 2-body 


BaBar 2-body 


BaBar Dalitz 




(140 fb" 1 ) 


(113 fb" 1 ) 


(192 fb" 1 ) 


A p7T 

A CP 


-0.16 ±0.10 ±0.02 


-0.114 ±0.062 ±0.027 


-0.088 ±0.049 ±0.013 




0.25 ± 0.17 t° f 6 


0.35 ±0.13 ±0.05 


0.34 ±0.11 ±0.05 




-0.28 ± 0.23 tom 


-0.13 ±0.18 ±0.04 


-0.10 ±0.14 ±0.04 




0.38 ± 0.18 


0.20 ±0.13 ±0.05 


0.15 ±0.11 ±0.03 




-0.30 ±0.24 ±0.09 


0.33 ±0.18 ±0.03 


0.22 ±0.15 ±0.03 


K- 


-0.02 ±0.16^2 
-0.53 ±0.29 to.ol 


-0.18 ±0.13 ± 0.05 
-0.52 tg;^ ± 0.07 


-0.21 ±0.11 ± 0.04 
-0.47 to.lt ± 0.06 




6 8 10 
At (ps) 



Figure 5. BaBar results for B°^p + n^ [24]: the 
At distributions of q = 1 tags (top), q = — 1 
tags (middle), and the resulting CP asymmetry 
(bottom). The smooth curves are projections of 
the unbinncd ML fit. 



4> 2 ; however, since the penguin contribution is 
unknown, additional information is needed. A 
recent theoretical model [28] uses SU (3) symme- 
try and the measured rates or limits for branch- 



AT^ttF, B° 



ing fractions of B° 
B ± ^K*°ir ± , and B ± ^p ± K°. S£/(3)-breaking 
effects are considered at tree level and accounted 
for via a factor f n /f K - The strong phase dif- 
ference between the two tree amplitudes (W — ► p 
and spectators p) is assumed to be small, as pre- 
dicted by factorization. The resulting central val- 
ues and errors for <f) 2 are: 102±19° for Belle values 
of C p7r , S p7r , AC p7r , AS pn ; 93 ± 17° for BaBar val- 
ues (113 fb- 1 ); and 95 ± 16° for Belle and BaBar 
values combined. 

The BaBar Dalitz plot analysis (192 ft)" 1 ) [26] 
allows one to directly fit for <f) 2 with little theoret- 
ical uncertainty from the penguin contribution. 
The result is <j> 2 = (113 (stat) ± 6(syst))°, 
consistent with the SU (3)-based results above. 



5. B c 



■p+p 



The decay B° — > p + p~ has two vector particles 
in the final state. If the p mesons are longitudi- 
nally polarized, i is even and CP = +1; but if 
they are transversely polarized, I can be even or 
odd and the final state is not a CP eigenstate. 

For longitudinal polarization, <j) 2 can be de- 
termined from the At distribution as done for 
B°^ir + 7r~. However, B° — > p + p~ has an advan- 
tage: the penguin contribution is expected to be 
small relative to the tree contribution [29] , which 
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reduces theoretical uncertainty on (j> 2 . Unfortu- 
nately B° — > p + p~ is more challenging experimen- 
tally: there are several backgrounds and also pos- 
sible nonresonant contributions. The method de- 
pends upon the p's being longitudinally polarized; 
otherwise a more involved angular analysis is nec- 
essary to determine (j> 2 [30]. Finally, the nonneg- 
ligible decay width of the p allows for I = 1 fi- 
nal states, which complicates extracting <p 2 via an 
isospin analysis [31]. 

The decay B° — > p + p~ has been observed by 
BaBar and the CP-violating parameters C pp and 
S measured with 81 fb _1 of data [7,8] and up- 
dated with 113 fb _1 of data [25]. A similar 
analysis is underway at Belle with 250 fb _1 of 
data. The final state consists of four pions, two 
charged and two neutral. In the case of multi- 
ple B° — >■ p + p~ candidates arising from multiple 



a, 



£,(» 



,00 



the p candidates 



m 7r o) is chosen, where i runs over 
From MC simulation, it is 
found that one or more pions from B° — > p + p~ 
are swapped with pions from the tag side 39% 
(16%) of the time for longitudinal (transverse) 
polarization. 

BaBar selects events with relatively loose cuts 
and does an unbinned ML fit to the At dis- 
tribution, including pdf's to account for back- 
grounds. Nonresonant contributions and inter- 
ference with decays yielding the same final state, 
e.g., B° — > a-yW , are estimated to be small and 
neglected. For 81 fb _1 of data, 24288 events are 
fit and a signal yield of 224 ± 29 is obtained. The 
fit includes a pdf for the angles 6 X 
i is the angle between the tt° from pf 
and the B° in the pf rest frame (i = 1, 2). This 
pdf has the form [32] 



and 9 2 , where 



d 2 T 



9 



T dcos9 1 dcost 



{n 



cos 2 6* n cos 2 # 



1-/ 



L ' sin 2 i 



a} (8) 



and determines f L , the fraction of longitudinally 
polarized decays. The fit results are C = 



-0.23 ±0.24 ±0.14 and S pp = 



-0.19 ±0.33 ±0.11 



(113 fb- 1 ), and f L = 0.99 ± 0.03 



±o :^ (si fb- 



G 

> 
pq 



20 
10 



r is 



G 

> 

pq 



a, 

(SI 

i 

'SI 

< 
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Figure 6. BaBar results for B° p + p~ [7]: the 
At distributions of q = 1 tags (top), q = — 1 
tags (middle), and the resulting CP asymmetry 
(bottom). The smooth curves are projections of 
the unbinned ML fit. 



systematic. 

It is fortunate that f L is close to unity; in 
this case the final state has CP — +1 and an 
angular analysis to determine <f> 2 is unnecessary. 
Fig. 6 shows the At distributions for q = ±1 
tagged events along with the resulting CP asym- 
metry. No CP violation is observed. Inputting 



the measured values for C„ 



and S pp into an 



isospin analysis that includes the branching frac- 
tions for B°^p+p- [7] and B+ -> p+ p° [33], and 
the upper limit for B(B° — * p°p°) [34], one obtains 
2 = (96 ± 10 (stat) ± 4(syst) ± ll theOT y)° [35]. 
The last error is due to the penguin contribution; 
it is significantly smaller than that for B — > mr 
(±36°), as expected. 

6. SUMMARY 

Time-dependent CP asymmetries in B° — > 



The first error listed is statistical and the second 



7T 7T 



and B () — > p + p decays are 
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measured and used to constrain the CKM an- 
gle <f>2- The B° -^tt + tt~ mode is experimentally 
clean but has the largest penguin contribution, 
which contributes theoretical uncertainty to <j> 2 . 
A model-independent constraint is 90° < (f> 2 < 
146° for \P/T\ < 0.45 (95% C.L.). An SU(3)- 
based model [18] indicates (f> 2 = (103 ± 17)° (and 
also that \P/T\ is large). Belle observes large CP 
violation in this mode while BaBar does not. 

The B° — > p^n^ mode is more complicated as 
there are more backgrounds than for £?°^7r + 7r~ 
and the final state is not a CP eigenstate. A 
model based upon SU (3) symmetry and using the 
measured branching fractions for B—tK*^ and 
B -> p ± K obtains <j> 2 = 95 ± 16° (Belle + BaBar 
quasi-two-body results combined). 

The B° — ► p + p~ mode has the smallest penguin 
contribution but suffers from additional back- 
grounds, possible nonresonant contributions, and 
a possible 1 = 1 component in the final state. 
Neglecting the latter two effects, BaBar mea- 
sures C pp and S pp for longitudinal polarization, 
which dominates the decay. Combining the mea- 
sured values with the branching fractions or lim- 
its for B" -> p+p- [7], B+ p+p° [33], and 
B" p°p" [34] gives (f> 2 = (96 ± 10 (stat) ± 
4 (syst) ± H thoory )° [35]. This value is similar to 
those obtained from measurements of B° -^ir + n~ 
and _B° — ► p ± n zf decays. 
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